Astronomy & Astrophysics manuscript no. sed 


February 2, 2008 


(DOI: will be inserted by hand later) 





Spectral energy distribution of the metagalactic ionizing 
radiation field from QSO absorption spectra* 

I. I. Agafonova 1 , M. Centurion 2 , S. A. Levshakov 1 , and P. Molaro 2 



in 
o 
o 

(N 
C 



> 

00 

m 

o 
in 
o 

^: 

Oh' 

6 

CO 



X 



1 Department of Theoretical Astrophysics, Ioffe Physico- Technical Institute, 194021 St. Petersburg, Russia 

2 Osservatorio Astronomico di Trieste, Via G. B. Tiepolo 11, 34131 Trieste, Italy 

Received 00 / Accepted 00 

Abstract. A computational procedure is presented to estimate the spectral shape of the ionizing background 
between 1 and 10 Ryd by analyzing optically thin absorption systems in the spectra of high redshift quasars. The 
procedure is based on the response surface methodology from the theory of experimental design. The shape of 
the recovered UV background at z ~ 3 shows a significant intensity decrease between 3 and 4 Ryd compared to 
the metagalactic spectrum of Haardt & Madau (1996). This decrease is interpreted as produced by Hell Gunn- 
Peterson effect. There are no features indicating a contribution from galaxies to the UV background which is 
therefore dominated by QSOs at z ~ 3. 
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1. Introduction 

The study of the spectral energy distribution (SED) in the 
radiation background resulting from flux emitted by all 
celestial sources is an important part of modern observa- 
tional cosmology. The current paradigm assumes that the 
metagalactic ionizing background is formed by radiation of 
QSOs and galaxies (stars) reprocessed by the intergalactic 
medium (IGM). The SED evolves with cosmic time due 
to different contribution of galactic and QSO radiation, 
and to varying IGM opacity caused by the hydrogen and 
helium reionization. 

Since the pioneering works of Chaffee et al. (1986) and 
Bergeron & Stasinska (1986), the approach to recover the 
shape of the ionizing radiation field from the measure- 
ments of the intervening metal absorbers in QSO spectra 
has been employed in numerous studies. In particular, the 
energy range 1 Ryd < E < 10 Ryd can be probed through 
the relative intensities of metal lines such as Sin-Siiv, 
C n-C IV, N in, N V, O VI. A common procedure consists of 
selecting a standard SED and checking whether it is con- 
sistent with measured column densities. This procedure 
can be appropriate for obtaining some general informa- 
tion about the SED, but its effectiveness and accuracy is 
rather low because of lack of search strategy. 
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In the present paper we describe a computational tech- 
nique exploiting the response surface methodology from 
the theory of experimental design which enables a directed 
search for the shape of the ionizing background. The pro- 
posed method reliably recovers the main features in the 
spectral shape of the metagalactic flux by the analysis of 
optically thin QSO absorption systems. Four absorption 
systems were selected to illustrate how this approach can 
be implemented in practice. 

Three of these systems have redshifts z ~ 3. This red- 
shift is of particular interest for studying the traces of 
still not completely ionized Hen (Reimers et al. 1997). In 
general, the presence of He II in the intergalactic medium 
affects the ionizing spectrum in the range E > 1 Ryd 
due to He n Lya, two-photon reemission and to He II con- 
tinuum absorption (Haardt & Madau 1996; Fardal et al. 
1998). The Hell Lya absorption troughs detected in spec- 
tra of five quasars (Vogel & Reimers 1995; Zheng et al. 
1998; Anderson et al. 1999; Heap et al. 2000; Smette et al. 
2002; Jakobsen et al. 2003 ; Zheng et al. 2004a) may indi- 
cate He II Gunn-Peterson effect (He n Lya absorption in a 
diffuse IGM) at z ~ 3. Furthermore, recent FUSE observa- 
tions of the He II Lya forest towards HE 2347-4342 com- 
bined with observations of the H I Lya forest at the VLT 
revealed also large scale variations in r\ = Hell/Hi which 
are still not well understood (Shull et al. 2004; Zheng et 
al. 2004b). One possible explanation is the presence of 
spatial fluctuations in the metagalactic ionizing field at 
E ^ 4 Ryd. To confirm this suggestion, additional obser- 
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vations at different rcdshifts and in different sightlines are 
needed. However, direct measurements of Hell Lya ab- 
sorption at 303.78 A are very problematic due to almost 
complete light blotting of distant QSOs by the intervening 
Lyman limit or damped Lya systems (LLSs and DLAs, 
respectively). In this context, the proposed approach to 
estimate the shape of the ionizing background from metal 
absorption systems is of great importance. 

The structure of the paper is as follows. The proce- 
dure to recover the shape of the underlying UV continuum 
and an example illustrating how to use it are described in 
Sect. 2. This section also contains a brief description of 
a computational method used to invert the observed line 
profiles - the Monte Carlo Inversion (MCI). The detailed 
analysis of absorption systems used for the SED estima- 
tions is given in Sect. 3. The results are discussed and 
summarized in Sect. 4. An example of an experimental 
design is given in the Appendix. 

2. Computational methods 

2.1. Shape of the ionizing radiation 

Here follows a description of how the SED of the ionizing 
radiation can be estimated from metal line profiles ob- 
served in optically thin absorption systems. The method 
is based on the response surface methodology used in ex- 
perimental design (see, e.g., Box, Hunter & Hunter 1978, 
Chapter 15). 

A basic UV spectrum is taken as an initial guess. The 
estimated column densities of different ions and a pho- 
toionization code (e.g., CLOUDY, Ferland 1997) are used 
to derive the ionization parameter, mctallicity and ele- 
ment abundance ratios. If the observed column densities 
are well reproduced with this trial spectrum, it implies 
that the initial guess was appropriate. If, however, the 
observed column densities cannot be reproduced or/and 
some other inconsistencies arise, then the assumed shape 
of the UV background is to be adjusted. 

To allow for quantitative estimations, the shape of 
the UV continuum must be parameterized. In general, 
the shape of the ionizing spectrum can be specified by 
a piecewise continuous function, e.g. by a set of power 
laws and/or exponents. In particular, in the range 1 Ryd 
< E < 10 Ryd relevant for the ions frequently observed 
in QSO spectra, the following variables (called 'factors' in 
the experimental design) can be used to describe the main 
features of the SED (Fig. 1): 

A first region between 1 Ryd and the Hei break is 
characterized by a slope /i (power law exponent) and by 
the position / 2 of the He I break. A second region can be 
defined between the He I and He n breaks (point A) with 
fz and /4 the slope and the position of the He n break, re- 
spectively. Since the energy depression can be relevant, we 
introduce here also f$ which defines the depth of the break 
[h = logC^WJA.)]) and fe - the slope between points A 
and B. f-j characterizes the energy spectrum in the re- 
gion B — C. To account for a possible bump around 3 Ryd 
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Fig. 1. Schematic picture of the z = 3 metagalactic (dashed 
line) and AGN-type (smooth line) ionizing backgrounds from 
Haardt & Madau (1996) and Mathews & Ferland (1987), re- 
spectively. The spectrum is normalized so that J„ {hv = 1 Ryd) 
= 1. The emission bump at 3 Ryd in the metagalactic spec- 
trum is caused by reemission of Hen Lya, Hen two-photon 
continuum emission and He II Balmer continuum emission from 
intergalactic clouds. The definition of factors {fi} is given in 
Sect. 2.1 



due to recombinations within the clumpy intergalactic gas 
(He II Lya emission, He n two-photon continuum emission 
and Hen Balmer continuum emission), we introduce fac- 
tors /s and /g which describe the amplitude of the bump 
and its width, respectively. 

The energy of the far UV cut off (point C) when taken 
above 100 Ryd does not affect the fractional ionizations of 
ions we are interested in. In all computations described in 
subsequent sections this energy and the slope after it are 
kept fixed at 128 Ryd and -1.5, respectively. Note that the 
positions of the He I and He n breaks can be shifted due 
to large scale motions and/or superposition of different 
spectra (e.g. stcllar+metagalactic), and hence, in general 
the values of ji and fi may not be equal to 1.8 Ryd and 
4 Ryd, respectively. Additional factors can be defined to 
describe some particular features of the SED as, for ex- 
ample, those caused by Hen Lya absorption in a diffuse 
IGM (see Sect. 3 below). 

In this notation, an ionizing spectrum is determined as 
a point {/i, . . . , fk} in the fc-dimensional factor space. The 
number and selection of factors depend on a particular 
spectral shape. For instance, the metagalactic spectrum 
at z — 3 calculated by Haardt & Madau (1996, hereafter 
HM) has coordinates {/i, / 4 , f 5 , fe, h, fs, M = { -1-4, 
0.55, -1.3, -13.0, -0.45, 0.46, 0.11} with /i standing for 
the entire slope between 1 Ryd and f^, whereas the point 
{/i,/2,/3,/4,/ 5 ,/6,/ 7 } = {'-0.5, 0.25, -1.0, 0.61, -2.5, 
—3.0,-0.7 } corresponds to the AGN-type spectrum of 
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Fig. 2. Velocity (upper panel) and gas density (lower panel) 
distributions along the line of sight used to produce the mock 
absorption spectra which are shown in Fig. 4 



Mathews & Ferland (1986, hereafter MF) (all energies are 
given in log 10 ). 

The next step of the shape adjustment procedure is to 
find a direction in the factor space which leads to the UV 
background with more appropriate characteristics. The 
starting ('null') point is represented by the initial spec- 
trum. A set of new trial spectra is produced by varying 
the factor values around the 'null' point in accordance 
with a special scheme called 'experimental (factorial) de- 
sign'. The experimental design is represented by a matrix 
with n (n > k) rows containing particular values of fac- 
tors (an example of the experimental design is given in 
Appendix) . 

To evaluate the fitness of each trial ionizing spectrum 
a numerical measure TZ (usually called 'response') is to 
be defined in such a way that bigger values of TZ corre- 
spond to increasing fitness. In general, the choice of TZ 
occurs rather heuristically and accounts for the informa- 
tion obtained with the initial UV spectrum and for any 
a priori information (like, e.g., allowable element abun- 
dance ratios). For illustration consider the following two 
examples. 

Assume that a set of silicon (Si II, Si III, Si iv) and car- 
bon (C II, C in, C iv) lines is observed in an absorption 
system, and that the initial UV spectrum overestimates 
the column densities of Si n and C iv and underestimates 
those of Si iv and Cn. Then the response can be written 
as 

- _ Cn Siiv 
U ~ CP X Sin' 

i.e., the search should go towards increasing the product 
of these ratios. 

Another example concerns an absorption system with 
the observed lines of Sim, Siiv, and Civ. The ratio 
Sim/Si IV determines the ionization parameter U. Using 




E, Ryd 

Fig. 3. A model ionizing background (smooth line) used to 
calculate the mock spectra shown in Fig. 4. As an initial ap- 
proximation for the restoring procedure an HM-type spectrum 
is taken (long-dashed line). The results of the first and second 
iterations are illustrated by the short-dashed and dotted lines, 
respectively 

the corresponding ion fractions Y, the ratio 
Si Siiv Tciv 

c = cPt^ 

can be calculated. Both observations and theoretical con- 
siderations give for this ratio a safe upper limit Si/C 
< 3(Si/C)©. If, for instance, the initial UV background 
delivers Si/C = 10(Si/C)O, then the search for a new SED 
can be governed by the ratio 



calculated for the ionization parameter U given by the 
observed column density ratio Sim/Siiv (note that the 
value of U depends on the shape of the trial spectrum). 
More examples are presented in the subsequent sections. 

The calculated responses {TZi}f =1 can be considered as 
points belonging to some surface in /c-dimcnsional factor 
space (referred to as 'response surface'). To determine the 
direction of steepest ascent the response surface should be 
described by some analytical function. A standard linear 
model is a reasonable first approximation: 

k 

K = Y,<*ih + P, (l) 

where {ai}f =1 are factor effects, /, is the scaled and cen- 
tered value of the ith factor, /j = (/j — /o,i)/si, «i is some 
suitable scale (range of variation) of the ith factor, and 
P = a a + Ylij=i a i,jfifj stands for the joint effect of the 
free term ao and of the nonlinear term that contains also 
the interaction of factors. 
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Fig. 4. Mock spectra of hydrogen and metal absorption lines (dots with 1 a error bars). The MCI solution corresponding to the 
initial ionizing spectrum is shown by the dotted curves. The smooth curves are the synthetic profiles obtained with recovered 
2nd iteration SED (dotted line in Fig. 3) 



The coefficients a, and (3 and their dispersions are cal- 
culated from n values of 1Z according to special formulas 
(depending on the type of experimental design used). The 
validity of the planar data model (1) can be checked by 
comparison of the estimation, (3' , with the value of IZq at 
the 'null point', i.e., when /j = ,/o.i- Statistical significance 
of the difference (IZo — /?') points to non-negligible non- 
linear effects. In this case the higher-order data model is 



to be used with corresponding higher-order experimental 
design. 

Given the factor effects {c^}, the maximization of 1Z 
is obtained by moving the factor values from the 'null' 
point in the direction normal to the response surface. This 
movement is performed stepwise until either the desired 
UV background is found (i.e. the one which reproduces 
self-consistently all observed column densities) or the va- 
lidity limit of the employed data model is reached. In the 
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Parameter 
(1) 


Model 
values 
(2) 


Estimated with 
initial SED 

(3) 


Estimated with 
1st iteration SED 

(4) 


Estimated with 
2nd iteration SED 

(5) 


u 


1.55E-2 


5.7E-2 


1.3E-2 


1.4E-2 


Nh, cm" 2 


1.2E19 


2.5E18 


8.4E18 


1.0E19 


<7 V , km s" 1 


20.0 


17.1 


22.2 


18.0 


y 


0.40 


0.56 


0.50 


0.48 


Zc 


8.0E-6 


2.7E-5 


1.0E-5 


9.0E-6 


Z S i 


1.0E-6 


3.8E-6 


1.6E-6 


1.2E-6 


iV(Hl), cm" 2 


3.4E15 


3.7E15 


3.4E15 


3.4E15 


iV(Cll), cm" 2 


2.5E12 


3.0E12 


2.9E12 


2.6E12 


AT(Cm), cm" 2 


8.1E13 


5.5E13 


6.9E13 


7.7E13 


AT(Civ), cm~ 2 


9.0E12 


9.9E12 


9.2E12 


8.6E12 


JV(Sill), cm" 2 


3.5E11 


4.2E11 


3.9E11 


3.5E11 


JV(Siin), cm" 2 


5.5E12 


5.2E12 


5.4E12 


5.3E12 


iV(Siiv), cm" 2 


1.9E12 


2.2E12 


1.8E12 


1.9E12 


Note: = Nx/Nh 



latter case the adjustment procedure should be repeated, 
this time with a newly obtained UV spectrum as a 'null' 
point. 

The initial set of factors may turn out to be redun- 
dant, i.e. some of the factor effects ai may be statistically 
insignificant. This means that the corresponding features 
in the spectral shape do not affect the absorption line pro- 
files included in the analysis. Such insignificant factors are 
removed from the data model (1) and are fixed at some 
appropriate level. The effects of the remaining m factors 
are then recalculated using the corresponding m-factorial 
design. 

The uncertainty of the recovered spectral shape re- 
quires some comments. The spectrum of the ionizing ra- 
diation, J„, defines the rates (s -1 ) of the photoionization 
processes, T i3 -, through the integral 

T H = j ^^-aii( v ) dv . ( 2 ) 

where Oij (v) is the photoionization cross section at fre- 
quency v for species i in ionization state j, and is the 
threshold frequency for photoionization. 

Thus, the SED is obtained by solving the integral equa- 
tion which is a classic ill-posed problem. Its solution de- 
pends crucially from the number of constraints included 
in the analysis. When applied to our problem, this means 
that the SED is best recovered if the corresponding ab- 
sorption system reveals many lines of different metals in 
different ionization stages. If only a few metal lines are 
observed, they can be used to recover the shape not in the 
entire range 1 Ryd < E < 10 Ryd, but in some narrower 
regions. For example, the lines of Cm, Civ and Ovi al- 
low to estimate the SED in the vicinity of the He II break 
(3 Ryd < E < 4.5 Ryd). 



2.2. The Monte Carlo Inversion procedure 

Absorption systems are analyzed by means of the 
Monte Carlo Inversion (MCI) procedure described in de- 
tail in Levshakov, Agafonova & Kegel (2000, hereafter 
LAK), and with modifications in Levshakov et al. (2002, 
2003a,b). Here we briefly outline the basics needed to un- 
derstand the results presented below in Sect. 3. 

The MCI is based on the assumption that all lines 
observed in the absorption system are formed in a contin- 
uous medium where the gas density, nn(x), and velocity, 
v(x), fluctuate from point to point giving rise to complex 
profiles (here x is the space coordinate along the line of 
sight). 

The MCI also assumes that within the absorber the 
metal abundances are constant, the gas is optically thin 
for the ionizing UV radiation, and the gas is in thermal 
and ionization equilibrium. The intensity and the spectral 
shape of the background ionizing radiation are treated as 
external parameters. 

The radial velocity v(x) and gas density nn(x) are 
considered as two continuous random functions which are 
represented by their sampled values at equally spaced in- 
tervals Ax. The computational procedure uses the adap- 
tive simulated annealing. The fractional ionizations of 
all elements included in the analysis are computed at 
every space coordinate x with the photoionization code 
CLOUDY. 

In the MCI procedure the following physical parame- 
ters are directly estimated: the mean ionization parame- 
ter Uq, the total hydrogen column density Vh, the line- 
of-sight velocity dispersion cr v , and density dispersion a y , 
of the bulk material [y = hh(x)/ho]! and the chemical 
abundances Z a of all elements involved in the analysis. 
With these parameters we can further calculate the col- 
umn densities V a for different species, and the mean ki- 
netic temperature T^ n . If the absolute intensity of the UV 
background is known, then the mean gas number density 
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Fig. 5. Hydrogen and metal absorption lines associated with the z a b s = 2.9171 system towards HE 0940-1050 (histograms). 
Synthetic profiles corresponding to the MCI solution with the z — 3 HM ionizing spectrum (smooth line in Fig. 6) are shown by 
the dotted curves, whereas those calculated with the final adjusted SED (short-dashed line in Fig. 6) are the smooth curves. The 
corresponding physical parameters are listed in Table 2, Col. 2. Bold horizontal lines mark pixels included in the optimization 
procedure. The zero radial velocity is fixed at z = 2.9171 



no, and the line-of-sight thickness L of the absorber can 
be evaluated as well. 

In general, the uncertainties of the fitting parameters 
U , Nn, and Z a are about 15%-20% (for data 

with S/N > 30) and the errors of the estimated column 
densities are less than 10%. However, in individual ab- 
sorption systems, the accuracy of recovered values can be 
lower due to different reasons such as partial blending of 
line profiles, saturation of profiles or absence of lines of 
subsequent ionic transitions. 



The procedure of spectral shape adjustment described 
in the preceding section is used in a modified form. From 
computational point of view, it is more convenient within 
the MCI to evaluate responses 1Z using ion fractions rather 
then column densities (at a given metallicity, the column 
density of an ion is proportional to its fraction). The pro- 
cedure stops when spectral shape ensuring x 2 < 1 for all 
lines observed in the system is found. 
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Fig. 6. SEDs recovered from the absorption line systems. The 
positions of ionization thresholds of different ions are indi- 
cated by tick marks. The definition of factors {/»} are given in 
Sect. 3.1. For details, see text 



2.3. Testing the SED recovering procedure 

In order to test the procedure of the SED adjustment we 
prepared a mock absorption line spectrum using the den- 
sity and velocity distributions plotted in Fig. 2, the UV 
ionizing background shown in Fig. 3 (solid line), and the 
model parameters listed in Table 1, Col. 2. The result- 
ing absorption lines after convolution of the intensities 
with a Gaussian-type point-spread function of FWHM = 
7 km s _1 and the addition of white noise with disper- 
sion 0.02 (S/N = 50) are shown in Fig. 4 by points with 
error bars. The chosen UV background can be produced 
by a mixture of stellar and metagalactic spectra (see, e.g., 
Giroux & Shull 1997). 

The procedure requires an input SED and as initial 
guess for the underlying UV background we took the 
metagalactic spectrum of HM shown as a long-dashed 
line in Fig. 3. This spectrum was defined by 7 factors 
{/i> A; /s> /bj hi fs, fo\ with the meanings described in 
Sect. 2.1. The only change is for /i which now stands for 
the entire slope between 1 Ryd and the Hell break at f&. 

The best MCI solution gives the physical parameters 
listed in Table 1, Col. 3. The corresponding synthetic ab- 
sorptions are shown in Fig. 4 by the dotted lines. As ex- 
pected, not all absorptions arc well reproduced, since the 
tried UV background differs from the one used to generate 
the mock lines. An underestimation of the C in and Si in 
intensities along with overestimation of C iv is well seen 
in Fig. 4. 

To calculate responses, the following expression was 
chosen: 

,5 , ( T C iii TgmA 
K = log — , 

V I CIV 1 Sill / 



with an additional constraint log(Tsnii/Tsiiv) <^ 0.4 
[i.e. the fractions of the corresponding ions are calcu- 
lated for the value of U determined by the condition 
log(T Si „i/T S iiv) > 0.4]. 

Factors were varied at 2 levels according to a 7-factorial 
saturated simplex design which gave 8 new trial UV 
spectra (see Appendix). These spectra were inserted into 
CLOUDY and with the obtained fractional ionizations 
the responses for every spectrum were evaluated. Then 
the factor effects were estimated, linearity of data model 
checked, and a new UV background produced by moving 
the factor values in the direction normal to the response 
surface. 

The new solution for the background is shown in Fig. 3 
by the short-dashed line, and the results of the MCI cal- 
culations are given in Table 1, Col. 4. Fitting of synthetic 
profiles to most absorption lines becomes much better, but 
C in, C iv and Si iv lines still have large x 2 values. Thus, 
the adjustment procedure was repeated once more, with 
the same set of factors, but different response 

t5 , /TCcrv T~siiv\ 
ft = log ~ ~ , 

V l CU 1 Sill / 

with constraints log(Tcm/Tciv) > 0.9 and 
log(T S iiii/T SiIV ) > 0.5. 

The resulting UV spectrum is shown in Fig. 3 by dots, 
and the corresponding synthetic profiles by the smooth 
lines in Fig. 4. The recovered parameters are listed in 
Col. 5, Table 1. Now all absorption lines are well repro- 
duced, and, thus, the inverse procedure can be considered 
completed. 

This test demonstrates both the effectiveness of the 
adjustment procedure and its limitations. Starting from 
some standard spectral shape, we were able to recover 
quite correctly all significant features of the underlying 
ionizing spectrum such as the depth of the Hen break 
and its shift to the lower energies. On the other hand, 
some fine features of the underlying spectrum (break at 
1.8 Ryd, small bump at 3 Ryd) cannot be reproduced with 
the adopted noise level. Higher S/N data would be more 
appropriate for this case. 

3. Application to observed absorption systems 

Here we analyze some QSO absorption systems using the 
procedure described in the preceding section. All compu- 
tations below were performed with laboratory wavelengths 
and oscillator strengths taken from Morton (2003). Solar 
abundances were taken from Asplund, Grevesse & Sauval 
(2005). 

3.1. Absorption system at z a ^ s = 2.9171 towards HE 
0940-1050 

This absorption system shown in Fig. 5 exhibits 
higher order Hi Lyman series lines and transitions of 
Cii/Ciii/Civ, and Si n/Si in/Si iv. 
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Fig. 7. Hydrogen and metal absorption lines associated with the z a bs = 2.9659 system towards Q 0347-3819 (histograms). 
Synthetic profiles corresponding to the MCI solution for the final adjusted SED (dot-short-dashed line in Fig. 6) are shown 
by the smooth curves. The corresponding physical parameters are listed in Table 2, Col. 3. Bold horizontal lines mark pixels 
included in the optimization procedure. The zero radial velocity is fixed at z — 2.9659 



The system was described in detail in Levshakov et al. 
(2003c) as an extremely low metallicity LLS with iV(Hi) 
= 3 x 10 17 cm~ 2 and Z = 0.001 Z Q . Levshakov et al. also 
noted that the observed line intensities were inconsistent 
with the HM metagalactic ionizing spectrum, and that 
better fitting of the observed profiles corresponded to the 
ionizing spectrum having a significantly enhanced Hen 
re-emission bump at 3 Ryd. This result was obtained by 
common 'trials and errors' method. Here we re-analyze the 
z a b s = 2.9171 system using the procedure of the directed 
search. 

With the HM ionizing spectrum at z — 3 as initial 
guess we obtained the synthetic profiles shown in Fig. 5 
by dotted curves. Strong overestimation of Sin A1260 A 
is clearly seen along with underestimation of C n A1334 A 
and Si iv A1393 A. Whether Cm A977 A line is underes- 
timated or not, cannot be decided unambiguously at this 
stage of investigation since the apparent intensity may be 
due to contamination by some Lya interloper. The HM 
spectrum was parameterized using the same 7 factors as 
in the example described in Sect. 2.3. The responses were 
calculated in the form 



K = log 



Ten T Si iv 



The design employed was the 7-factorial saturated simplex 
design (see Appendix). The first-iteration ionizing spec- 
trum revealed a sharp He n break shifted to 3 Ryd and the 
He II re-emission bump significantly increased as compared 



to the initial approximation (Fig. 6, dotted line). This 
spectrum provided better fitting to the observed profiles 
(e.g., all carbon lines including Cm were now described 
with x 2 ~ 1), but the profiles of Sin and Siiv remained, 
correspondingly, over- and underestimated. 

Next iteration of the SED adjustment was performed 
with the same set of factors, design and response aug- 
mented with constraints log(Tcin/Tciv) < 1 and 
log(Tsiin/Tsiiv) > 0.1 to enable self-consistent descrip- 
tion of the C in and Si ill lines. However, this step failed 
to produce any improvement. Attempts to include addi- 
tional factors describing the Hei break (see Fig. 1) or to 
employ other forms of the responses failed as well. 

The solution was found after performing model cal- 
culations (using CLOUDY) of radiation transmission 
through a plane-parallel absorbing cloud. For an HM-type 
ionizing spectrum, r\ = iV(He n)/A^(H i) > 50 and the ab- 
sorber with A^(Hi) = 3 x 10 17 cm -2 , marginally opti- 
cally thick in hydrogen continuum, is certainly opaque in 
Hen. The radiation in the Hen continuum is effectively 
absorbed by such cloud and a part of this radiation is 
reemitted as the Hen A304 A Lya and two-photon ree- 
mission, and He n Balmer continuum emission. As a result, 
the incident HM spectrum being transmitted through the 
absorber with A^(Hi) = 3 x 10 17 cm~ 2 reveals a sharp 
and deep break at 4 Ryd and a strong emission line at 3 
Ryd of about twice the incident intensity. Thus, the in- 
creased intensity at E < 3 Ryd in the ionizing spectrum 
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Fig. 8. Same as Fig. 7 but for the 2 a bs = 2.9375 system towards HE 0940-1050. The MCI synthetic profiles (smooth curves) are 
calculated with the final adjusted SED shown by the dot-long dashed line in Fig. 6. The corresponding physical parameters are 
listed in Table 2, Col. 4. The zero radial velocity is fixed at z = 2.9375 



recovered from the z a bs = 2.92 absorber is, probably, due 
to auto-emission of the cloud. 

However, in the recovered spectrum a strong intensity 
depression starts just above 3 Ryd, and this result has 
been reproduced in all trials with different responses and 
different sets of factors. At the same time, the model cal- 
culations show that the spectral range 3 Ryd < E < 4 
Ryd remains unaffected by the processes inside the cloud, 
i.e., if the intensity depression in this range is not present 
in the incident spectrum, then it does not appear in the 
transmitted one. This means that the depression at E > 3 
Ryd in the ionizing spectrum is produced by the processes 
outside the cloud. We can assume that this is a manifes- 
tation of He II Lya A304 A absorption arising in both 
smoothly distributed IGM (Hen Gunn-Peterson effect) 
and discrete Lya forest clouds (line blanketing) . He II Lya 
absorption in the diffuse IGM gas results in the absorp- 
tion trough blueward of the resonant wavelength 304 A 
(3 Ryd). This trough together with the subsequent ab- 
sorption in the Hen continuum (E > 4 Ryd) produce a 
winding structure in the spectral shape at E > 3 Ryd. 
As a first approximation, this structure can be described 
by a straight step (see Fig. 6). 

After this consideration, a new set of factors was de- 
fined to account for this step-like form (Fig. 6): f\ - the 
slope between 1 and 3 Ryd; f 2 - the depth of the Hen 
absorption through; / 3 - the energy of the He n ionization 
break; fi - the depth of the Hen ionization break; f$ - 
the slope of the Hen ionization break; fe,fr - the height 
and width of the Hen re-emission bump; / 8 - the slope 
after the He n ionization break. In the following analysis, 
factors and /§ turned out to have low effects and were 
set to 4 Ryd and -0.45, respectively. 



The final UV background is shown by the short-dashed 
line in Fig. 6. The synthetic line profiles calculated with 
this SED are plotted by the smooth solid lines in Fig. 5 
and the corresponding physical parameters are listed in 
Table 2, Col. 1. The predicted column density of Hell for 
this absorber is 2.8 x 10 19 cm" 2 giving r\ — 87. 

In the framework of our approach, the restored spec- 
trum may be considered as some average spectrum. This 
spectrum reveals features stemming from both the inter- 
galactic incident radiation (depression between 3 and 4 
Ryd in Fig. 6) and local processes in the cloud itself (en- 
hanced re-emission of Hen Lya, deep break at 4 Ryd). 
To reconstruct the true metagalactic spectrum, absorbers 
optically thin in He II are needed, i.e. those with AT (Hi) 
< 10 15 cm" 2 . 

3.2. Absorption system at = 2.9659 towards Q 
0347-3819 

This system has been described in detail in Lcvshakov et 
al. (2003b). Since then Q 0347-3819 was re-observed at the 
VLT/UVES with a full wavelength coverage and longer 
exposure time. This allowed to obtain a high-quality spec- 
trum with resolution of ~ 6 km s _1 and S/N = 50 — 100. 
Here we repeat the analysis for the new data. 

The hydrogen and metal absorption lines of the z a b s = 
2.9659 system are shown in Fig. 7 (histograms). Neither 
Cn A1334 A (undetected), nor Sin A1260 A (blended) 
can be used. This makes the presence of C in A977 A line 
crucial for the SED estimation since Sim, Siiv and Civ 
lines can be fitted with a very broad range of ionizing 
spectra (constrained only by the condition [Si/C] < 0.5). 
Although we did not find any metal line candidate for 
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Fig. 9. Same as Fig. 7 but for the z a b s = 1.9426 system towards J 2233-606. The MCI synthetic profiles (smooth curves) are 
calculated with the SEDs shown by the dotted and dashed lines in Fig. 10. The corresponding physical parameters are listed in 
Table 2, Col. 5. The zero radial velocity is fixed at z = 1.9426 



blending with C III, we cannot exclude blending with a 
Lya forest absorption. Thus, the analysis below should be 
taken with caution and it is valid only if the absorption 
at the position of C in is entirely due to this ion. 

Calculations with the HM ionizing background showed 
that this spectrum significantly underestimates C III. Our 
trials to reconstruct the underlying continuum shape were 
conducted with two sets of factors, with and without ac- 
counting for the Hen Lya absorption. In both cases it 
was possible to obtain the appropriate ionizing spectra. 
However, the spectrum estimated without Hen Lya ab- 
sorption shows significantly enhanced emission at 3 Ryd 



and at the same time is harder at E > 4 Ryd compared to 
the initial HM spectrum. Since the absorber under study 
has iV(Hi) = 3.4 x 10 14 cm~ 2 and is optically thin in Hell, 
this type of ionizing spectrum is obviously an artifact and 
should be rejected. 

The recovered spectrum with the intensity depression 
at 3 Ryd< E < 4 Ryd is shown in Fig. 6 by the point- 
dashed line with the corresponding synthetic profiles plot- 
ted by the smooth lines in Fig. 7 (physical parameters are 
given in Table 2, Col. 3). Since the depth of the Hen Lya 
trough (factor /2) affects quite strongly the ion fractions, 
its value is estimated with accuracy better than 0.1 dex. 
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The depth of the Hen continuum absorption (factor f±) 
has an uncertainty of about 0.15 dex. The predicted Hen 
column density is 2.1 x 10 16 cm~ 2 (rj = 62) and, hence, 
this cloud is optically thin in Hen, t^"* ~ 0.03. Thus, 
the recovered background coincides in this case with the 
infalling radiation spectrum and is not modified by local 
effects. From the depth of the Hell Lya trough we can 
estimate the opacity of the intergalactic He 11 Lya absorp- 
tion as THell = 2.1 ± 0.2. 

3.3. Absorption system at z a b s — 2.9375 towards HE 
0940-1050 

The absorption system at z a b s = 2.9375 towards HE 0940- 
1050 reveals metal lines of C in A977 A, C iv AA1548, 1550 
A and Ovi A1037 A (Ovi A1031 A is blended, as well 
as Nv AA1238, 1242 A) (Fig. 8) . We cannot exclude that 
the observed Ovi A1037 A and Cm lines are contami- 
nated by some hydrogen absorption. However, no metal 
line candidates for blending were found. As in the case of 
the z a bs = 2.9659 system described in the preceding sub- 
section, the analysis below is valid only if the intensities 
arc entirely due to absorption of the corresponding ions. 

Two lines of one element along with one line of another 
can be described with a very broad range of ionizing spec- 
tra. However, different spectra will deliver different abun- 
dance ratios [O/C] and this can be used to constrain the 
allowable SEDs. Measurements of [O/C] in Galactic and 
extragalactic H 11 regions and in metal poor halo stars in- 
dicate that a safe upper bound for [O/C] is 0.5 (Henry, 
Edmunds, & Koppen 2000; Akcrman et al. 2004). 

Calculations with the HM spectrum produce for this 
absorption system [O/C] > 1 suggesting that the shape 
of the input ionizing background could be inadequate. To 
meet [O/C] < 0.5, the ionizing spectrum should be either 
significantly harder (by 0.5 dex) at E > 4 Ryd or should 
have a step-like depression at E > 3 Ryd. The second 
option is preferable since this step is also present in the 
spectrum recovered from the absorber at z a b s = 2.9171 
along the same line of sight (Sect. 3.1). The step at 3 
Ryd < E < 4 Ryd strongly affects the Cm/Civ ratio 
leading to the result that the observed ratio is obtained 
at a higher ionization parameter as compared to the HM 
spectrum and, hence, the corresponding Ovi fraction is 
larger. 

Using the HM spectrum as an initial approximation 
and varying the depth of the He 11 Lya absorption trough 
(i.e. one-factor experiment) we adjusted the spectral shape 
in a way that enables the description of the observed lines 
with [O/C] < 0.5. The obtained spectral shape is shown 
in Fig. 6 (long-dashed line), with the synthetic profiles 
plotted by the smooth lines in Fig. 7 and the physical 
parameters given in Table 2, Col. 4. The uncertainty in 
the step depth is of 0.1 dex (given the energy above 4 Ryd 
at the level of the HM spectrum). The predicted column 
density of the once ionized helium A^(Hell) = 8.2 x 10 16 




Fig. 10. SEDs recovered from the absorption line system at 
z = 1.9426 towards J 2233-606. For details, see text 



cm 2 translates to 77 = 146 which is almost 2.5 times 
higher as in the z a b s = 2.9659 system. 

3.4. Absorption system at z.^ s = 1.9426 towards J 
2233-606 

This system has been described by Prochaska & Buries 
(1999), D'Odorico & Petitjean (2001) and Levshakov et 
al. (2002). A large wavelength coverage provided by com- 
bining the VLT/UVES and HST/STIS spectra, allows to 
identify many higher order hydrogen lines and numerous 
metal transitions (Fig. 9). Unfortunately, the important 
Cn A1334 A line is blended with a Ly-a forest absorp- 
tion, and C 11 A1036 A is very noisy, but this system is 
nevertheless worth being analyzed due to the rare occa- 
sion of simultaneous presence of low- (Si 11, Mg 11) and high 
ionization lines like N v. 

Calculations with the z = 2 HM ionizing spectrum 
(Fig. 10, solid line) describe most of the observed line pro- 
files except Nm A989 A and Nv AA1238, 1242 A which 
come over- and underestimated, respectively. The SED 
was adjusted in two experiments with factor sets both 
accounting for the Gunn-Peterson (GP) absorption and 
without it, and with the response maximizing the product 
of the ratios T N v/T Nm , T C ii/T C iv, and T S iiv/T Sim . 
All trials showed that to reproduce the nitrogen lines 
the ionizing background should be significantly harder at 
E > 4 Ryd than the initial HM spectrum. Unfortunately, 
because of a low quality of the C 11 AA1334, 1036 A lines, 
it turned out to be impossible to distinguish between the 
ionizing spectra either with or without the GP absorption. 
In fact, both recovered SEDs (shown by the dotted and 
dashed lines in Fig. 10) give identical fitting to the data. 
The synthetic profiles are shown by the smooth lines in 
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Table 2. Physical parameters of the metal absorbers derived by the MCI procedure (limits are given at the 1 a level) 



Parameter" 


z abs = 2.9171 
HE 0940-1050 
( z ) 


z abs = 2.9659 

Q 0347-3819 


z abs = 2.9375 

HE 0940-1050 
(V 


z abs = 1.9426 
J 2233-606 


U 


2.46E-2 


2.26E-2 


1.87E-1 


1.97E-2 


JVh, cm -2 


6.00E20 


7.22E17 


3.19E19 


3.20E19 


<r v , km s" 1 


28.4 


14.2 


35.2 


33.5 


°y 


0.75 


0.54 


0.70 


0.75 


Zc 


3.8E-7 


1.02E-4 


3.0E-6 


5.3E-5 


Zn 


< 1-3E-8 


9.0E-6 




4.9E-6 


Zo 






1.2E-5 


< 3.0E-4 


Zyig 








1.2E-5 


Zai 


< 3.8E-9 






8.0E-7 


Z S i 


4.9E-8 


2.9E-5 




9.9E-6 


[Zc] 


-2.8 


-0.38 


-1.90 


-0.66 


[Zn] 


<-3.7 


-0.82 




-1.1 


[Zo] 






-1.57 


<-0.2 


[Zm s ] 








-0.44 


[Zai] 


^ -2-8 






-0.52 


[Zsi] 


-2.8 


-0.05 




-0.51 


TV(Hl), cm -2 


3.2E17 


3.4E14 


5.6E14 


2.5E16 


iV(Cll), cm -2 


5.8E12 


<9.0E11 




4.0E13 


iV(Mgll), cm -2 








5.3E12 


iV(Sill), cm -2 


8.9E11 






5.6E12 


iV(Cni), cm" 2 


2. OEM 


4.8E13 


5.7E12 


1.0E15 


iV(NlIl), cm" 2 


< 7.8E12 


4.3E12 




9.8E13 


iV(Allll), cm" 2 


< 2.2E11 






1.5E12 


TV(Silll), cm" 2 


1.2E13 


3.8E12 




7.9E13 


iV(Clv), cm" 2 


2.4E13 


1.8E13 


9.2E12 


3.9E14 


AT(Silv), cm" 2 


7.6E12 


6.0E12 




5.8E13 


iV(Nv), cm" 2 








9.2E12 


AT(Ovi), cm" 2 






4.2E13 


< 1.4E14 


(T),K 


2.8E4 


1.3E4 


3.7E4 


1.8E4 


Af(Heii) 6 , cm" 2 


2.8E19 


2.1E16 


8.2E16 


1.1E18 


v b 


87 


62 


146 


44 



a Z x = Nx/Nh; [Z x ] = log(Nx/Nu) - lo g (Nx/N n ) Q . 
''Predicted values. 



Fig. 9. However, an upper limit can be set for a putative 
GP Hell absorption at z = 2: th ii < 1-8. 

The predicted column density of Hen of 1.1 x 10 18 
cm" 2 (r] = 44) gives the optical depth in the He II contin- 
uum of 1.5. In principle, taken at face value such optical 
depth can soften noticeably the incident ionizing contin- 
uum at E > 4 Ryd. However, the observed lines of Nm 
and Nv clearly require a hard spectrum in this energy 
range. This contradiction can be explained in two ways: 
cither the Nv lines arise in the external regions of the 
absorbing cloud where the incident ionizing continuum is 
not yet distorted by the He n continuum absorption, or the 
density variations along the line of sight make the effec- 
tive optical depth smaller. In both cases we can conclude 
that the metagalactic ionizing spectrum at z ~ 2 is harder 
at E > 4 Ryd as compared to the spectrum predicted by 
HM. 



4. Conclusions 

We have proposed a method to estimate the shape of the 
ionizing continuum from metal lines observed in the inter- 
vening absorption systems. The implementation includes 
the following steps: 

- parameterization of the spectral shape by means of a 
set of factors based on physical processes relevant for 
the IGM; 

- choice of some quantitative measure (response) to es- 
timate the fitness of a trial shape; 

- generation of trial shapes by randomization of the fac- 
tor values according to some experimental design; 

- evaluation of the factor effects; 

- estimation of the optimal shape by moving in factor 
space towards the ascending fitness. 

The result depends on the number of metal lines involved 
in the analysis: in general, the more lines of different ionic 
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transitions of different elements are detected in an absorp- 
tion system the higher is the accuracy of the recovered 
spectral shape. In some cases additional information con- 
cerning, for instance, element abundance ratios can signif- 
icantly tighten the factor values, in particular when only 
a few metal lines are available. 

Although the main objective of this work is to illus- 
trate how the proposed approach can be used in practice, 
some physical results are worth mentioning as well. They 
are as follows: 

1. The metagalactic ionizing spectrum at redshift z ~ 3 
has breaks at 3 and 4 Ryd and is quite hard at E > 4 
Ryd, at least as hard as a model metagalactic spectrum 
of Haardt & Madau (1996). 

2. The intensity decrease between 3 and 4 Ryd is prob- 
ably produced mostly by He n Lya absorption in the 
intergalactic diffuse gas since line-blanketing from dis- 
crete Lya forest clouds cannot account for significant 
He II opacity with such hard ionizing spectrum (Fardal, 
Giroux & Shull 1998; Zheng, Davidsen & Kriss 1998). 
Thus, the intensity depression at 3 Ryd < E < 4 Ryd 
in the spectrum of the metagalactic ionizing radiation 
may be an imprint of a true He n Gunn-Peterson effect. 

3. The He n opacity estimated from the depth of the He II 
Lya absorption trough is th ii ~ 2. This should be 
considered as an average optical depth since in our 
approach the Hen Lya absorption is approximated 
by a straight step. As known from observations, Hen 
opacity at z ~ 3 fluctuates revealing both absorption 
troughs and opacity gaps. For reference, th c ii ~ 4 is 
measured in absorption troughs at 2.77 < z < 2.87 
in the spectrum of HE 2347-4342 (Zheng et al. 2004) 
along with THoii < 0.5 in the opacity gaps at z a bs = 
2.817 and 2.866 (Reimers et al. 1997). 

4. A putative galactic contribution to the ionizing back- 
ground spectrum (e.g. Pettini et al. 2001; Ciardi, 
Bianchi & Ferrara 2002; Fujita et al. 2003), if present, 
would significantly soften the spectrum at E > 4 Ryd. 
We do not find any traces of soft component in the 
recovered spectra: even in the absorption systems at 
z = 2.96 and z = 1.94, where high metallicity supposes 
their kinship with galaxies, the observed line intensi- 
ties point to a hard ionizing background. Thus, we can 
conclude that the ionizing background at z ~ 3 and 
z ~ 2 is dominated by QSOs. 
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Appendix A: An example of experimental design 
for 7 factors 

The factors are tested at two levels: fi — fos ± en, where 
o i is a suitable scale. Usually Oi ~ (0.1 — 0.2)/ ,i, in order 
to maintain the linear dependence of the response TZ on 
the factors. The normalized and centered factor values at 
the upper and lower levels are, respectively, +1 and — 1. 
Then a 7-factorial saturated simplex design can be defined 
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by the following symmetric and orthogonal matrix {fij} 
(Nalimov 1971): 





fl 


h 


h 


h 


h 


h 


A 


ii 


+1 


-1 


-l 


-l 


+i 


+i 


+i 


-i 




+1 


+1 


-l 


-l 


-l 


-l 


+i 


+i 


n 2 


+1 


-1 


+i 


-l 


-l 


+i 


-l 


+i 


n 3 


+1 


+1 


+i 


-l 


+i 


-l 


-l 


-i 


n 4 


+1 


-1 


-l 


+i 


+i 


-l 


-l 


+i 


K 5 


+1 


+1 


-l 


+i 


-l 


+i 


-l 


-i 


■Re 


+1 


-1 


+i 


+i 


-l 


-l 


+i 


-i 


n 7 


+1 


+1 


+i 


+i 


+i 


+i 


+i 


+i 


n 8 



The first column is added to estimate the free term (3 in the 
model (1). Such experimental design is called a 'two level 
saturated plan' since here all degrees of freedom N = 8 
are used to estimate 8 regression coefficient ({ai}J =1 and 

1 N 



and 



3 = 1 



1 N 



Other plans with different numbers of levels and treat- 
ments (rows in the matrix) can be used as well. 

After the regression coefficients for the model (1) are 
found, new factor values can be calculated from 



fi 



fo,i + AK^fii , 



where \ii = ai/ ^Jj2j=i a ] arc direction cosines and A1Z 
is the increment of the response. The value of A1Z is re- 
stricted by the condition x 2 ~ 1 for every absorption line 
(or its portion) involved in the optimization procedure. 



